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Abstract 

This paper describes flavour independent searches for hadronically decaying 
neutral Higgs bosons in the data collected by the DELPHI experiment at LEP, 
at centre-of-mass energies between 189 and 209 GeV. The collected data-set 
corresponds to an integrated luminosity of around 610 pb _1 . The e + e~— > hA 
and e + e~— > hZ processes are considered, with direct Higgs boson decays into 
hadrons. No evidence for Higgs boson production is found, and cross-section 
limits are set as a function of the Higgs boson masses. No explicit assump- 
tions are made on the underlying physics beyond the Standard Model, allowing 
interpretation of the data in a large class of models. 
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1 Introduction 

As long as bosonic decay channels are kinematically closed, the Standard Model (SM) 
Higgs boson decays preferentially into the heaviest accessible fermion pair. This leads to 
the dominance of Higgs boson decays to b quarks in most of the mass range accessible 
at LEP. In extensions of the SM, however, the Higgs boson couplings to b quarks might 
well be suppressed. This can for example occur in the Minimal Supersymmetric Standard 
Model (MSSM) [1], and in general Two Higgs Doublet Models (2HDM) [2]. Suppressed 
couplings to b-quarks are also possible in models with composite Higgs and gauge bosons, 
where the dominant Higgs boson decay channel is gluonic [3]. The SM Higgs boson search 
at LEP strongly relies on the identification of b quarks or r leptons to separate a possible 
signal from background processes, and has reduced sensitivity to the final states predicted 
by such alternative models. 

This paper presents searches that only assume direct hadronic decays of the Higgs 
bosons. The scope of these analyses is thus more general, and encompasses the models 
referred to above. However, rather than stating the results in explicit theoretical frame- 
works, model independent bounds on production cross-sections are derived, allowing to 
confront a wide class of models with the data. 

Two Higgs boson production processes are studied in this paper: the Higgsstrahlung 
process e + e~— > hZ, and Higgs boson pair production e + e~— > hA, which is present in 
all extensions of the Standard Model Higgs sector beyond one doublet. Although the 
notation is familiar from CP conserving models, the CP properties of h and A are not 
necessarily specified here. 

Flavour independent searches for hadronic Higgs boson decays have been published 
by the other LEP collaborations [4-6] . 

The paper is organized as follows. The common framework and search strategy is 
described in Section 2, and is followed by an overview of the experimental setup and 
the data sets in Section 3. The search for hadronic decays of hA pairs is described in 
Section 4. Section 5 describes the search for hZ production, where the Higgs boson 
is accompanied by jets, missing energy, electrons, or muons from the Z boson decay. 
Section 6 combines the various analyses, and gives flavour independent results on the hZ 
and hA cross-sections, as a function of the Higgs boson masses. Section 7 concludes the 
paper. 

2 General features 

The analyses presented in this paper rely essentially on kinematic and event shape 
information. Without the background rejection allowed by b quark or r lepton iden- 
tification, the signal contamination by quark pair production and/or vector boson pair 
production can be important, and dijet mass information must be used to discriminate 
signal from background. The way in which this information is exploited depends on the 
search channel, and is detailed in the corresponding sections. 

The search results independence from the hadronic flavour of the Higgs bosons decay 
products is obtained by considering only the decay mode with lowest sensitivity in each 
search channel. The precise procedure varies again from one channel to the other, and 
details are given for each analysis in the corresponding description. 

All search results presented in this work are interpreted using a modified frequentist 
technique based on the extended likelihood ratio [7]. For a given experiment, the test 
statistic Q is defined as the likelihood ratio of the signal+background hypothesis (s+b) 



2 



to the background hypothesis (b), computed from the number of observed and expected 
events in both hypotheses. Measured discriminating variables can be used to assign a 
signal to background ratio to individual events. Probability density functions (PDFs) for 
Q in the b and s+b hypotheses are built using Monte Carlo sampling of the expected 
background and signal rates, and of the optional discriminating variable distributions. 
The confidence levels CL b and CL s+b are defined as the integrals of the b and s+b PDFs 
for Q between — oo and the observed value Q bs- 

A signal hypothesis to which there is no sensitivity might be uncorrectly excluded by 
CL s+b in case a downward fluctuation of the data occurs. To prevent this, we adopt the 
conservative approach of defining the confidence level in the signal hypothesis, CL S , as 
the ratio CL s+b /CL b . 1-CL S measures the confidence with which the signal hypothesis 
can be rejected, and will be larger than 0.95 for an exclusion confidence of 95%. 

When no signal is found, upper bounds on the production cross-section times the 
branching fraction into hadrons are extracted for both production processes as a func- 
tion of the Higgs boson masses. These bounds are expressed in terms of reference 
cross-sections, defined as follows. The production rate of any final state resulting from 
e + e~— > hZ can be expressed in terms of the SM hZ cross-section, a™ [8]. The reference 
cross-section for e + e - ^ hA, <rj^{, is calculated with HZHA [9], assuming the 2HDM and 
the absence of any mixing in the Higgs sector. The result depends only on electroweak 
parameters and the h and A Higgs boson masses, and remains valid for models with any 
number of Higgs doublets or singlets. Suppression terms coming from branching fractions 
and possible suppression of the Higgs boson couplings to the Z (such factors are hereafter 
denoted by BR and R, respectively) all factorize, and are hidden in a generic (model 
independent) suppression factor : 



Ohz^z+hadrons = &lz x R ^z x BR(h -> hadrons) 

— °hZ x u hZ) 

°"hA-»hadrons = °hA x -^ hA x B R(h — > hadrons) x BR(A — > hadrons) 



Our results will be expressed in terms of C bZ and C^ A . 

Although the traditional CP-conserving Higgs boson nomenclature is used throughout 
the paper, part of the results can also be extended to the context of CP-violating Higgs 
sectors. For hA production, the parity and charge conjugation properties of the Higgs 
bosons do not play any role in the decay kinematics. The results on C^ A are therefore 
universally valid for the pair production of any two Higgs bosons (i.e. hihj — > hadrons), 
regardless of their CP properties. On the contrary, our results on e + e~^> hZ assume 
standard quantum numbers for the Higgs boson. Non-standard Higgs boson parity prop- 
erties would imply different polarization of the associated Z boson, thereby affecting the 
angular distributions of the bosons and the selection efficiencies. Our results for this 
channel should therefore be used with care. 

If kinematically open and allowed by the Higgs boson quantum numbers, cascade 
decays like h— >AA may become dominant. The efficiency of the analyses developed 
below to these final states has not been evaluated. 

The study also applies to non-Higgs scalar particles. The cross-sections and analyses 
presented here however assume that the produced scalars have negligible width (less than 
about 1 GeV/c 2 ). 
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3 Detector, data samples and simulation 

A detailed description of the DELPHI detector and its performance during the LEP1 
data taking period can be found in [10]. For LEP2, the vertex detector was upgraded [11], 
and scintillator counters were added in regions without electromagnetic calorimetry (at 
40 and 90 degrees from the beam axis), improving the detection of events with radiative 
photons in these regions. The detector was operated in nominal conditions until Septem- 
ber 2000, when one of the twelve sectors of DELPHI'S main tracking device underwent an 
irremediable failure. A special algorithm, based on complementary subdetectors, allowed 
to limit the tracking efficiency loss in this region. 

The data samples used in this paper were collected during the last 3 years of LEP 
operation (1998 to 2000), and are clustered around seven centre-of-mass energy values, 
as shown in Table 1. The total used luminosity amounts to 610.4 pb _1 . The hZ analysis 
uses the data from all listed sets, while the hA search uses all but the two sets with lowest 
luminosity. 



year 


1998 


1999 


2000 


x/i (GeV) 
C (pb- 1 ) 


188.6 


191.6 


195.5 199.5 


201.6 


205.0 206.5 


158.0 


25.9 


76.9 84.3 


41.1 


82.0 142.2 



Table 1: Centre-of-mass energy and integrated luminosity of the data collected by the 
DELPHI detector in the years 1998-2000. 

The simulation samples used for this study have been produced at the centre-of-mass 
energies given in Table 1; details are given below. In all cases, samples have been pro- 
cessed that comply with the experimental situation before and after the detector failure 
mentioned above. 

For the hZ channel, two-fermion background events are simulated using the KK2F [12] 
generator, while WPHACT [13] is used for four-fermion final states (see [14] for the specific 
DELPHI implementation). For the hA analysis, prepared during an earlier stage of 
this work, the two-fermion and four-fermion SM backgrounds are simulated using the 
PYTHIA [15] and EXCALIBUR [16] generators respectively. Backgrounds from two-photon 
events are generated using PYTHIA. All background samples have a size equivalent to 
20-50 times the luminosity recorded at the corresponding centre-of-mass energy. 

The HZHA generator is used to generate hZ signal samples for 40<m h <120 GeV/c 2 , 
in steps of 2.5 GeV/c 2 . Z boson decays into quarks, neutrinos, electrons and muons are 
simulated, and the Higgs boson is decaying to either a pair of gluons, or a pair of light 
quarks (s quarks in practice); at several mass points, samples with Higgs boson decays into 
heavy quarks are used to verify the flavour-independence of the results. For the neutrino 
channel, signal samples are also generated every 5 GeV/c 2 for 10<mh<40 GeV/c 2 , and at 
mh=4 GeV/c 2 ; heavy quark samples are simulated at all mass points (when kinematically 
allowed) in this case. All hZ signal samples contain 5000 events. 

The hA process is also simulated using the HZHA generator, in the mass range 
16<m h +m A <180 GeV/c 2 . The simulation starts at m h ,m A =4 GeV/c 2 ; steps of 5 GeV/c 2 
are used when 10<rrih,mA<30 GeV/c 2 ; 10 GeV/c 2 steps are used otherwise. Higgs boson 
decays into gluon or light quark pairs are simulated over the whole mass range; samples 
with heavy quarks are simulated at several mass points for cross-checks. All hA signal 
samples contain 2000 events. 
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For low Higgs boson masses, the decay final state is determined by non-perturbative 
mechanisms (direct coupling to hadronic resonances). There is no robust model of this 
regime, and the efficiencies of the analyses described below are difficult to predict. We 
therefore limit our search to Higgs boson masses above 4 GeV/c 2 , where the PYTHIA 
hadronization model is expected to be valid. 

4 Search for hA production 

The flavour-independent hA search is designed to cover a large part of the accessible 
h and A mass range, and is based on general kinematic features such as event shapes and 
dijet mass information. 

The higher multiplicity of gluon jets compared to quark jets results in a higher selection 
efficiency for hA — > gluons, but also in a worse mass resolution for this signal. To 
minimize biases that may arise from these competing effects, the selection efficiencies are 
determined using samples of hA — > quarks, while the mass reconstruction is done with 
samples of hA — > gluons. 

4.1 Analysis streams 

A first preselection is applied to all events, requiring at least 20 charged particles, a 
total reconstructed energy greater than 60% of the centre-of-mass energy, and an effective 
centre-of-mass energy after initial state radiation greater than 150 GeV (the method used 
to estimate the effective centre-of-mass energy is described in [17]). The efficiencies of 
the multiplicity cut are typically 98% for the hA — > quarks signal samples, and 100% for 
the hA — > gluons samples. 

In the rest of the analysis, three different sets of cuts are designed to obtain sufficient 
signal acceptance over a large domain in the (mh,mA)-plane. The analysis streams are 
described below. 

4.1.1 Four-jet stream 

Close to the kinematic limit and when both Higgs bosons have comparable masses, a 
four-jet topology is expected. To analyse this topology, events are clustered into four jets 
with the DURHAM algorithm [18]. All jets are required to have an invariant mass larger 
than 2 GeV/c 2 , and to contain at least two charged particles. Events are retained in this 
stream if their thrust value is below 0.85 and if the product of the smallest jet energy 
and inter-jet angle (called E m i n a m i n ) is greater than 10 GeV-rad. Dijet invariant mass 
information is used to reject events compatible with WW production as in [19], requiring 
the corresponding probability, Pww> to be less than 0.01. This proves helpful not only 
when m h = m A ~ 80 GeV/c 2 , but also for other signal masses, where WW production 
contributes to the expected background through wrong jet pairings. 

4.1.2 Three-jet stream 

When the h and A mass difference increases, the events tend to contain three jets, since 
the jets resulting from the decay of the lighter Higgs boson cannot always be resolved. 
The same behaviour is observed if h and A both have low masses, because of the larger 
boost in this case. To analyse this topology, only events with thrust values between 0.70 
and 0.92 are kept. Events are then clustered into four jets, and events compatible with 
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WW production are rejected as in the four-jet stream. The remaining events are clustered 
into three jets. As before, all jets are required to have a mass larger than 2 GeV/c 2 and 
to contain at least two charged particles. 

4.1.3 High-thrust stream 

Finally, when both Higgs bosons are sufficiently light, signal events tend to become 
cigar-like. Events are selected in this analysis stream if they have a thrust value larger 
than 0.92, and are clustered into three jets, each jet having to satisfy the same quality 
criteria as above. As the dominant background comes from two-quark processes, the 
kinematic compatibility with W-pair events is not tested. 

4.1.4 Comparison between data and simulation 

Distributions of thrust, E min ai m i n and — log 10 Pww are shown in Figure 1 for example 
signals representative of the three analysis streams. Their performance and comple- 
mentarity are illustrated in Table 2. As discussed above, these efficiencies are obtained 
assuming h and A decays into light quark pairs. From comparisons with a number of 
samples with Higgs boson decays into b quarks, the quark flavour-dependence of the ef- 
ficiencies is found smaller than 2%. The efficiencies obtained with gluonic Higgs boson 
decays are higher by a factor up to 1.2 in the three-jet and four-jet analyses, and by a 
factor up to 2 in the high-thrust analysis. A numerical comparison between data and 
expected background is shown in Table 3. 



Signal 


tour-jet /(y \ 
e hA I/O ) 


three-jet /q/\ 
e hA { /Oj 


high-thrust / 07 n 
e hA I/O ) 


Four-jet events 


~ 70 


~ 45 


0-5 


Three-jet events 


~ 50 


~ 70 


0-5 


High-thrust events 


0-5 


0-5 


10-15 



Table 2: hA search. Efficiencies of the three analysis streams on three classes of hA signal 
events, considered as four-jet events when m^jfn^ > 60 GeV/c 2 ; as high-thrust events 
when m h ,mA < 30 GeV/c 2 ; and as three-jet events in the remaining cases. 



centre-of-mass energy 


189 GeV 


196 GeV 


200 GeV 


>204.5 GeV 


Four-jet stream : 










expected background 


433.5 


221.3 


259.0 


634.5 


observed events 


459 


248 


232 


642 


Three-jet stream : 










expected background 


1593.3 


750.6 


797.0 


1894.3 


observed events 


1585 


736 


772 


1824 


High-thrust stream : 










expected background 


1384.9 


642.5 


654.3 


1516.3 


observed events 


1331 


612 


607 


1450 



Table 3: hA search. Numerical comparison between data and background simulation, 
after all selections. The year 2000 data are merged in the last column. The statistical 
uncertainty on the background estimates is 2%. 
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Figure 1: hA search. Comparison between data and simulation for the variables used in 
the three selection streams, at the preselection level. Thrust, E m i n a: m i n , and -log 10 (Pww) 
are shown from top to bottom. The signal normalization assumes C^ A =1. The dis- 
criminating power is indicated by corresponding distributions from representative signal 
samples. In particular, the thrust distributions justify the use of three complementary 
selection streams. 
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4.2 Final discriminant variable 

In all three analysis streams, a four-constraint fit is performed on the selected events, 
requiring total energy and momentum conservation. The compatibility of each event with 
a given mass hypothesis (toi^toa) is estimated using the following quantity : 

2 ( m\ ec — \ 2 / m r 2 ec — m\ 

X {m h ,m A ) = — - rpr + 



In the above expression, rri[ ec is the mass of a given dijet, is the mass of the opposite 
dijet in the four-jet stream, and of the opposite jet in the other streams, and are the 

corresponding errors as obtained from the kinematic fit. For every event, the jet pairing 
that minimizes x 2 ( m h,mA) is used to compute the discriminant. 

The distributions of this variable are computed, for data and background, on 
a 1 GeV/c 2 x 1 GeV/c 2 grid in the (mh, mA)-plane. 

The signal distributions are first determined at the simulated mass points. To con- 
struct the signal distributions at mass configurations that are not simulated, the his- 
tograms are linearly interpolated between the three closest simulated points. The validity 
of this procedure is discussed in the next section. 

These distributions are used for the statistical evaluation of the compatibility of the 
data with the simulation following the procedure described in Section 2. Figure 2 illus- 
trates this procedure for three representative example signals. 



4.3 Systematic uncertainties 

Systematic uncertainties affecting this analysis are related to imperfect simulation of 
the detector and of SM processes, to residual flavour-dependence of the selections, and 
to the signal interpolation procedure. All effects are discussed below. 

The numerical differences at the preselection level between the observed data and the 
expectation from SM processes are taken as a first contribution to the uncertainty on 
the background prediction. For every subsequent selection variable, the difference of the 
average values of the corresponding distributions in data and background is computed; 
the selection cut is then varied by this amount. This procedure leads to an estimated 
systematic uncertainty on the backgrounds of 5% in the four-jet and three-jet streams, 
and 10% in the high-thrust stream, at each centre-of-mass energy. The higher uncertainty 
in the high-thrust stream is mainly due to the observed numerical discrepancy after the 
preselection. 

As stated earlier, the flavour- dependence of the selection efficiencies is estimated by 
comparing the results obtained for gluon decays, light quark decays, and b-quark decays. 
The preselection efficiency is slightly worse (by at most 2%) in the case of light quark 
decays, and is applied to the gluon decays as well. The remaining selections have a 
quark flavour dependence smaller than 2%. This number is taken as a contribution to 
the systematic uncertainty on the signal rate. 

Possible uncertainties from the interpolation procedure are illustrated in Fig. 3, which 
compares a simulated distribution with the result of an interpolation between neighbour- 
ing points. The shapes are well reproduced. To quantify more precisely what accuracy 
can be expected from this procedure, the previous exercise is repeated on a large number 
of simulated mass configurations. In each case, the bin contents of the simulated and 
interpolated histograms are compared. The distribution of the ratio of true and inter- 
polated bin contents has a spread of 5%, which is taken as a second contribution to the 
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Figure 2: hA search. Comparison between data and simulation for the distributions 
of x^ijn^mx). The discrimination obtained for (m^mA) = (70,60) in the four-jet 
stream, for {m^mx) = (100,20) in the three-jet stream, and for (m^rriA) = (30,20) 
in the high-thrust stream are shown from top to bottom. In all distributions, overflows 
have been included in the last bin. All the data used for this search are included (cf. 
Section 3). The signal normalization assumes C^ A =1. 
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systematic uncertainty on the signal rate. Note that this is a conservative estimate as 
the real interpolation is over smaller intervals than the 20 GeV/c 2 used here. 
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Figure 3: hA search. Left: distribution of x 2 (H0, 40), as obtained from an interpola- 
tion between the nearest simulated samples at (mh,mA)=(H0,30) and (110,50) GeV/c 2 . 
It is compared with the true (dashed line) distribution from a simulated sample at 
(mh,mA) = (H0,40). Right: distribution of x 2 (H0,30) for Higgs bosons decaying into 
either qq (solid line) or gg (dashed line), reflecting the worse mass resolution obtained 
for a four-gluon final state. In all distributions, overflows have been included in the last 
bin. 

Figure 3 also shows the discriminant variable distribution obtained for quarks and 
gluons, at (m h ,mA) = (110,30). As expected, the worse mass resolution on gluon jets 
translates into a somewhat reduced discriminating power; this is true for all mass points. 
The present analysis, using a selection efficiency determined on light-quark signal samples 
and a kinematic analysis calibrated on gluon samples, thus has no bias towards a given 
flavour that would invalidate the results when considering other Higgs boson hadronic 
decays. 

5 Search for hZ production 

This section describes a flavour-independent search for the e + e - -^ hZ process. The 
analyses used in this section are adapted from those used for the DELPHI e + e~— > ZZ 
and e + e~^ Z7* cross-section measurements [19,20], or from the previous section. The 
analyses are briefly described in the following; we refer to the original publications for 
further details. 

5.1 Final states with jets 

The hZ search in the fully hadronic final state is described below. Two mass re- 
gions are analysed separately, namely m h <40 GeV/c 2 using a three-jet analysis, and 
m h >40 GeV/c 2 using a four-jet analysis. 
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5.1.1 Low mass search (ra h < 40 GeV/c 2 ) 

The region 4 GeV/c 2 <mh<40 GeV/c 2 is analysed using the hA three-jet stream, de- 
scribed in section 4.1. Efficiencies and discriminant variable distributions are determined 
from interpolations between existing hA signals, obtaining distributions of x 2 ( m h,^z), 
where m z is the Z boson mass and m h is ranging from 4 to 40 GeV/c 2 . 

A few hZ samples have been simulated (at =4, 20 and 40 GeV/c 2 ) to verify the 
validity of this procedure. Despite the different Higgs boson angular distributions and 
the Z boson width, the efficiencies and x 2 ( m h,^z) distributions proved accurate within 
the systematic errors estimated in the hA three-jet channel. 

5.1.2 High mass search (ra h > 40 GeV/c 2 ) 

The selection used to analyze the four-jet channel in the e + e~— > TIL cross-section 
measurement [19] is adapted to test the hypothesis of Higgs boson production. After the 
selection of high multiplicity events with at least four DURHAM jets, the only remaining 
backgrounds are from ZZ and WW four-quark final states, and four jet events arising 
from qq production with hard final state gluon radiation. Event shape variables and jet 
kinematics are used to further enhance the signal. 

A constrained fit is performed on the reconstructed jets, requiring energy and momen- 
tum conservation between the initial and final state. The fitted jet momenta are then 
used to reconstruct the dijet masses for all possible jet pairings. The uncertainties on 
the fitted momenta and their correlation matrix are taken into account in that proce- 
dure. The dijet masses and their uncertainties are used to compute the probability that 
each event is compatible with WW, ZZ, qq, or hZ production (in this case, the event 
probability depends on the signal mass hypothesis), as in [19]. 

Finally, for each jet, the probability that it originates from a b quark is computed, 
using the methods described in [21]. Adding this information in the characterization of Z 
boson dijets allows better sensitivity to the fraction of the signal where the Z boson decays 
into b quarks. This information is not used on the Higgs boson side, thus preserving 
flavour- independence . 

The method employed to combine the various mass, event shape and b-tagging prob- 
abilities using likelihood ratio products is described in [19]. The resulting discriminant 
variable is called Pmggs- 

Table 4 gives a numerical comparison between data and simulation, for a few Higgs 
boson masses. Examples of ln(P H igg S ) distributions are shown in Figure 4. At each Higgs 
boson mass hypothesis, the Pmggs distributions in data, background and signal are used 
for the statistical interpretation of the results (see Section 2). For every mass hypothesis, 
the signal distribution is obtained from an interpolation between the two closest simulated 
signal samples. 

The systematic uncertainties are dominated by the uncertainty on the four-jet rate 
in qq events [22], conservatively taken to be 10% over the full Higgs boson mass range. 
Other contributions, evaluated in detail in [19], are small and can be neglected. 

5.2 Final states with jets and missing energy 

The event topology in the missing energy channel changes significantly as a function 
of the Higgs boson mass. For low enough Higgs boson masses a mono-jet final state is 
expected, while at high mass the jets are well separated, resulting in a two jet final state. 
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Figure 4: hZ search, final state with jets. Comparison between data and simulation for 
the distributions of ln(PHi ggs )- Higgs boson mass hypotheses of 75, 90 and 105 GeV/c 2 
are shown from top to bottom. All data from 189 to 209 GeV are included. The signal 
normalization assumes C\ z =\. 
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m h (GeV/c 2 ) 



SM (no hZ) observed 



€ hZ (%) hZ(h^qq) 



40 
50 
60 
70 
80 
90 
100 
110 



635.0 659 

522.8 532 

1784.1 1824 

5457.3 5476 

1681.7 1764 

957.4 970 

368.5 372 
87.6 75 



31.3 154.6 

30.1 132.3 
51.9 197.1 
78.6 249.3 

54.2 136.0 
54.6 96.2 
47.0 41.1 
33.9 9.4 



Table 4: hZ search, final state with jets. Numerical comparison between background 
simulation, data, and a few example signals, after an illustrative cut on Pmggs that max- 
imises the product of signal efficiency and purity. All data from 189 to 209 GeV are 
included. The signal normalization assumes C^ z —1. The statistical uncertainty on the 
signal efficiency and on the background is less than 1%. 

For both mass domains (with a transition set at m h =40 GeV/c 2 ), dedicated analyses are 
described below. 

5.2.1 Low mass search (m h < 40 GeV/c 2 ) 

The analysis used to measure the Z7* — > uuqq cross-sections [20] is adopted without 
major modifications to test the hypothesis of low mass Higgs boson production in the 
missing energy channel. After a preselection reducing the 77, Bhabha and Z7 back- 
grounds, three selections are designed, focusing on different expected topologies as a 
function of the mass of the Higgs boson. 

The first selection is optimized to probe very low Higgs boson masses, using an explicit 
cut on the visible mass of the event, which is required to be below 6 GeV/c 2 . The second 
selection exploits the large energy imbalance of huu final states: events are split into two 
hemispheres according to the plane perpendicular to the direction of the thrust axis, and 
are required to have one of the hemispheres containing at least 99% of the total visible 
energy in the event. The third selection, which is less efficient at very low masses because 
of an explicit cut on the charged particle multiplicity, is mainly based on a topological 
requirement: events are forced into a two jet configuration and an upper cut on the 
opening angle of the two jets is set at 78 degrees. All three selections use the information 
from the veto counters, by rejecting events with hits in veto counters far away from energy 
depositions in calorimeters or reconstructed tracks. The three analyses are combined on 
an event-by-event basis, by selecting events that pass at least one cut. 

After all selections, the most important remaining background processes are Z7* and 
Wei/. Higgs boson signal efficiencies reach up to 65%, and drop to 40% for masses below 
5 GeV/c 2 and above 40 GeV/c 2 . The efficiencies are found to be independent, within 
their statistical uncertainties, of the quark or gluon nature of the Higgs boson decay. The 
reconstructed visible mass distribution is used as discriminant variable for the statistical 
evaluation of the compatibility of the data with the expectation for different Higgs boson 
mass hypotheses (see Section 2). For every mass hypothesis, the signal distribution is 
obtained from an interpolation between the two closest simulated signal samples. Higgs 
boson decays to gluons are assumed, given the worse mass resolution expected in this 
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case. The reconstructed mass spectra and Higgs boson signal efficiencies are illustrated 
in Figure 5. 
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Figure 5: hZ search, final state with jets and missing energy, low mass search. Left: com- 
parison between data and simulation for the hadronic invariant mass distribution, after 
the mono-jet uuqq event selection. All data from 189 to 209 GeV are included. Right: se- 
lection efficiency as a function of the Higgs boson mass, at y/s = 206.5 GeV. 



The systematic uncertainty on the predicted background rate is 9%, obtained by vary- 
ing each cut by an amount equal to the difference of the corresponding variable's average 
values in data and in simulation, at the preselection level. The uncertainty also in- 
cludes a 5% contribution from the statistical uncertainty on the simulated background 
sample [20]). The uncertainty on the signal selection efficiency is 5%. 

5.2.2 High mass search (m h > 40 GeV/c 2 ) 

The selection described below closely follows the search for invisibly decaying Higgs 
bosons produced together with Z bosons decaying into hadrons [23]. Backgrounds origi- 
nating from qq(7) and 77 processes are rejected by applying the same selections as in [23]. 
Events with more than two identified leptons are then rejected. After this step, all events 
are clustered into both two-jet and three-jet configurations, using the DURHAM cluster- 
ing algorithm. 

To obtain a good performance over the mass range 40<mh<115 GeV/c 2 , three mass 
windows are defined (40<m h <70 GeV/c 2 , 70<m h <90 GeV/c 2 , and 90<m h <115 GeV/c 2 , 
subsequently referred to as first, second, and third mass range, respectively), and each 
is treated with a dedicated analysis. The discrimination is obtained in each case using 
a two-step Iterative Discriminant Analysis (IDA) [24]. After the first iteration, a cut on 
the IDA output is made to remove part of the background while retaining 90% of the 
signal. The IDA is then trained again on the events passing this selection. At each step, 
the IDA is trained simultaneously on ss, cc and bb flavours, resulting in a comparable 
performance for all flavours. 

Twelve variables compose the IDA discriminant for the first mass range. These are 
the scalar and vector sums of the transverse momenta of all reconstructed particles, the 
total visible energy, the event thrust, the energy of the least energetic jet in the three- 
jet configuration, the difference between the Fox- Wolfram moments H2 and H4 [25], 
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the energy of the most isolated particle, the complement of the two jets opening angle 
(acollinearity), the angle between the two planes formed by the jet axes and the beam 
axis (acoplanarity), the total missing mass, and the highest transverse momentum of a 
particle with respect to its jet in the two-jet configuration. The IDA was trained with 
signal samples ranging from mh=40 to 67.5 GeV/c 2 . 

The analysis in the second mass range uses the same twelve variables as in the low 
mass analysis, except for the energy of the least energetic jet in the three-jet configuration, 
which is replaced by the event b-tag probability as defined in [21]. The IDA was trained 
with signal samples ranging from mh=70 to 87.5 GeV/c 2 . 

In the third mass range, the IDA is composed of the same variables as above, but 
cuts removing the tails in the distributions of the input variables are applied prior to 
the training, as is done in [23] . For every centre-of-mass energy, the IDA was trained on 
signal samples ranging from mh=90 GeV/c 2 to the kinematic limit. 

Numerical comparisons between data and simulation are given in Table 5, after the 
cut on the first IDA output. 

The measured jet momenta are rescaled so that the missing mass equals mz, and the 
hadronic mass is computed using the rescaled jets. For each signal hypothesis, a window 
around the nominal Higgs boson mass is defined as a final selection, to maximize the 
product of signal efficiency and purity. 

The second step discriminant variable, called IDA2, is used in the statistical inter- 
pretation of the analyses. For every mass hypothesis, the signal distribution is obtained 
from an interpolation between the two closest simulated signal samples. The IDA2 dis- 
tributions are shown in Figure 6, for events inside the final mass window. This figure also 
shows that a better separation between background and signal is obtained for gluonic 
Higgs boson decays. 

The estimation of the systematic uncertainties also follows what was done in the search 
for invisible Higgs boson decays [23]. The dominant contribution originates from the 
description of the event shape variables, which affects the estimation of the qq(7) back- 
ground. Calibration data taken at centre-of-mass energies close to the Z pole are used 
to evaluate these systematics using the method described in [19,23]. Compared to these, 
the uncertainties on the background level due to the choice of jet clustering algorithm 
and the b-tagging procedure are small and can be neglected. 

The uncertainties on the background vary from year to year due to different detector 
operating conditions, and are between 5.7% and 12.6% for the first mass range, between 
6.2% and 10.2% for the second mass range, and between 4.0% and 8.2% for the third 
mass range. The corresponding systematic uncertainty on the signal efficiency is between 
1 and 3%, the largest uncertainty affecting the first mass range. 

5.3 Final states with jets and isolated leptons 

The qqe + e~ and qq/i + /i~ channels are analysed in the same way as in the e + e~— ■> ZZ 
and Z7* cross-section measurements. This search selects events with two identified leptons 
of the same flavour (electrons or muons) and of opposite charge, and at least two jets 
reconstructed with the JADE algorithm [26] (with y min = 0.01). Compared to [19,20], 
cuts on the invariant mass of the leptonic and hadronic systems (used to enrich the ZZ 
and Z7* samples) are not applied here, and the particle identification criteria are slightly 
modified in order to maximize the Higgs boson signal to background ratio. 

The reconstructed mass of the hadronic system is used as discriminant in the statis- 
tical interpretation of the analyses. For every mass hypothesis, the signal distribution 
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Figure 6: hZ search, final state with jets and missing energy, high mass search. Com- 
parison between data and simulation for the distributions of IDA2, after all selections. 
The distributions correspond to Higgs boson mass hypotheses of 75, 90 and 105 GeV/c 2 
from top to bottom, respectively. All data from 189 to 209 GeV are included. The signal 
normalization assumes C\ z =\. 
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Table 5: hZ search, final state with jets and missing energy, high mass search. Numerical 
comparison between background simulation, data, and a few example signals, just before 
the second IDA iteration. All data from 189 to 209 GeV are included. The signal 
normalization assumes C^ z =l. The statistical uncertainty on the signal efficiency and on 
the background is less than 1%. 

is obtained from an interpolation between the two closest simulated signal samples, and 
a window is defined around the nominal Higgs boson mass to maximize the signal to 
background ratio. 

Distributions of the reconstructed hadronic mass are shown in Figure 7. Tables 6 and 7 
display the signal efficiencies, and a numerical comparison between data and simulation. 

The main systematic uncertainties affecting the results are the uncertainty in the 
lepton identification and the signal selection efficiency. Combined with the statistical un- 
certainty from the finite size of the simulated sample, these result in a relative systematic 
uncertainty in selecting qqe + e~ and qq/i + yU~ events of 5%. The uncertainty on the small 
residual background rate is dominated by the MC statistics and amounts to 15%. 
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Figure 7: hZ search, final state with jets and isolated leptons. Comparison between data 
and simulation for the reconstructed hadronic mass distribution, when the mass of the 
opposite lepton pair is within 15 GeV/c 2 of the Z mass. All data from 189 to 209 GeV 
are included. 
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m h (GeV/c 2 ) 


SM(no hZ) 


observed 


e hZ (%) hZ(h^qq) 


40 


1.6 


3 


59.8 


14.2 


50 


1.1 


2 


61.8 


12.9 


60 


1.5 


3 


63.3 


11.4 


70 


2.3 


4 


63.4 


9.5 


80 


17.6 


18 


66.1 


7.9 


90 


23.3 


18 


69.5 


5.8 


100 


16.3 


11 


68.2 


2.8 


110 


1.9 


3 


64.2 


0.8 



Table 6: hZ search, final state with jets and isolated electrons. Numerical comparison 
between background simulation, data, and a few example signals, after the mass cuts 
that maximize the separating power for each Higgs boson mass hypothesis. All data from 
189 to 209 GeV are included. The signal normalization assumes — 1- 



m h (GeV/c 2 ) 


SM(no hZ) 


observed 


e hZ (%) hZ(h^qq) 


40 


1.5 


3 


77.1 


18.1 


50 


2.3 


4 


77.5 


16.1 


60 


2.0 


3 


78.5 


14.1 


70 


3.3 


3 


76.8 


11.5 


80 


20.0 


16 


78.0 


9.3 


90 


25.6 


20 


80.6 


6.7 


100 


17.8 


18 


81.1 


3.3 


110 


3.1 


2 


76.8 


0.9 



Table 7: hZ search, final state with jets and isolated muons. Numerical comparison 
between data, background simulation, and a few example signals, after the mass cuts 
that maximize the separating power for each Higgs boson mass hypothesis. All data from 
189 to 209 GeV are included. The signal normalization assumes C^ z —1. 
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6 Results 

The hypothesis of Higgs boson production is tested in this section, using the results 
from the analyses described above. The definition of the reference cross-sections, the sta- 
tistical procedure, and the scope and limitations of the results are discussed in Section 2. 
The results are represented graphically in this paper, but explicit exclusion values can be 
obtained from the Collaboration on request. 

6.1 Excluded cross-sections for hA production 

For each (mh,mA) hypothesis, the x 2 ( m h, %) distributions for background simulation, 
data and signal (see Figure 2) are used as input to the statistical analysis. 

The three analysis streams are combined as follows. The high-thrust stream is in- 
dependent from the three-jet and four-jet streams, and is statistically combined with 
each of them in turn. The resulting combinations partly overlap, and cannot be added 
as independent channels. The ambiguity is solved at each mass point by choosing the 
combination that provides the strongest expected exclusion. 

The mass domain is scanned from 771^,171 A >4 GeV/c 2 to the kinematic limit, in mass 
steps of 1 GeV/c 2 in both directions. At each mass point, the confidence level in the 
background hypothesis is computed, and the value of C^ A that is excluded with a con- 
fidence level of exactly 95% is determined. No evidence for signal is observed. Figure 8 
displays regions excluded in the (mh,mA)-plane, for a few fixed values of C^ A . For com- 
parable Higgs boson masses and maximal production cross-section (C 2 A =1), the combi- 
nations with mii+rriA below ~140 GeV/c 2 are excluded. When one of the Higgs bosons 
is very light, the mass of the other boson is constrained to be either below 4 GeV/c 2 
(the hA search threshold), or above 108 GeV/c 2 , still when C 2 A =1. The three-jet and 
high-thrust analysis streams allow a significant portion of the mass plane to be excluded 
even when C^ A is as small as 0.25. 

In most regions of the mass domain, the present search has high background and 
relatively low signal (see Figure 2). The exclusions are thus very sensitive to fluctuations 
in the data, which explains the irregularities in the exclusion regions. 

The above results are valid in the case of direct Higgs boson decays into hadrons. 
When the intermediate h— >AA decay mode opens (the corresponding kinematic domain 
is delimited by the dashed lines in Figure 8), the six-parton final state may become 
dominant. The sensitivity of the present analyses to this topology has not been evaluated. 

As stated in Section 2, the hA exclusions are valid independently of the parity prop- 
erties of h and A, and in particular may be interpreted in models with additional Higgs 
doublets or singlets, or in the context of a CP-violating Higgs sector. 

6.2 Excluded cross-sections for hZ production 

For each m h hypothesis, the discriminant variables provided by the analysis channels 
described in Section 5 are used as input to the statistical analysis. For masses below 
40 GeV/c 2 , these are the x 2 (^h,^z) distributions from the hadronic channel (like in 
Figure 2), and the mass spectra from the low mass missing energy channel (see Figure 5). 
Above 40 GeV/c 2 , the hadronic channel provides distributions of ln(PHi ggs (^h)) (Fig- 
ure 4), the high mass missing energy search provides distributions of the second step IDA 
output, IDA2 (Figure 6), and the channel with jets and leptons provides the hadronic 
mass distributions. 
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Figure 8: Exclusion regions in the (m h ,m A )-plane, for 0.5 < C^ A < 1, 0.25 < C£ A < 0.5, 
and C^ A < 0.25 from light to dark grey. The dashed lines indicate the borders where the 
h — > AA decay opens kinematically. The exclusions are given at the 95% CL. 
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The Higgs boson mass domain is scanned from mh>4 GeV/c 2 to the kinematic limit. 
At each tested mass, the confidence level in the background hypothesis is computed, and 
the value of Cj 2 z excluded with a confidence level of exactly 95% is determined. The scan 
is done in steps of 1 GeV/c 2 over most of the range, and in steps of 100 MeV/c 2 when the 
excluded value of C^ z approaches 1. Cross-section limits are computed for Higgs boson 
decays into either gluon or quark pairs, and the flavour-independent result at each mass 
is given by the weakest expected exclusion. 

Figure 9 shows the excluded values C\ z as a function of the Higgs bosons mass. Ob- 
served and expected limits agree well over a wide range of masses. The largest discrepancy 
occurs for a Higgs boson mass close to 30 GeV/c 2 where an excess of 2.65 a is observed 
(the excess is around 2.5 a when the result is averaged over the mass resolution at that 
particular Higgs boson mass hypothesis). The excess is shared by two channels: the 
missing energy channel contributes around 2 a, and the low-mass fully hadronic channel 
contributes around 1.5 a. The confidence level in the background hypothesis is shown 
as a function of mj, in Figure 10. The confidence level in the background reaches a few 
10~ 3 around m h =30 GeV/c 2 , but it should be remembered that given the large number 
of tested mass hypotheses, the probability that such an excess occurs at at least one mass 
is much larger. 

Except for signal hypotheses around mj, = 107 GeV/c 2 , there is a deficit of data in 
the hypothesis range 100 < m h < 115 GeV/c 2 . The discrepancy is at the level of 1 to 
1.5 a, as can be seen in Figure 10. The exclusions are thus stronger than expected in this 
region. Table 8 gives the mass limits for hadronic Higgs boson decays when C^ z =l, and 
for comparison the results achieved in the SM Higgs boson searches [27]. 





expected limit 


observed limit 




(GeV/c 2 ) 


(GeV/c 2 ) 


DELPHI (h^qq) 


108.0 


110.6 


DELPHI (h-*gg) 


109.2 


111.0 


DELPHI (SM decay) 


113.3 


114.1 



Table 8: The expected and observed 95% CL lower limits (in GeV/c 2 ) on the mass of a 
hadronically decaying Higgs boson, assuming its cross-section is identical to that in the 
SM. The results of the SM Higgs boson search are also given for comparison. 



7 Summary 

Searches for hadronic Higgs boson decays have been performed, exploiting the LEP2 
collision data recorded by DELPHI. No evidence of the presence of a signal was found in 
either the e + e~— »hZ or the e + e~— >hA process, in a mass range extending from 4 GeV/c 2 
to the kinematic limit. 

For the hZ process, a cross-section larger than 15% of the SM value (i.e C^ z >0.15) is 
excluded from m h =4 to 60 GeV/c 2 . Under the assumption of a production cross-section 
equal to that of the Standard Model, m h <110.6 GeV/c 2 is excluded at the 95% CL, while 
a limit of 108 GeV/c 2 is expected on average from background-only experiments. The 
cross-section limit is weaker than expected around mh=30 GeV/c 2 . 

A large part of the mass domain available for the hA process was explored as well. 
In the case of full production strength and 100% branching fraction into hadrons (i.e. 
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Figure 9: Excluded values of C^ z as a function of the Higgs boson mass, assuming direct 
hadronic Higgs boson decays. The solid curve is the observed exclusion, the median 
exclusion expected in the absence of signal is shown by the dashed curve, and the bands 
represent the 68% and 95% confidence intervals. 
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Figure 10: Confidence level in the background hypothesis as a function of the Higgs 
boson mass. The solid curve is the observed confidence, the median expectation from 
background-only experiments is shown by the dashed line, and the bands represent the 
68% and 95% confidence intervals. 
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Cl A =l), the excluded region extends roughly up to m h ,m A =108 GeV/c 2 for the heavier 
Higgs boson, when the lighter one has a mass below 10 GeV/c 2 . When both masses 
are equal, m h ,m A <70 GeV/c 2 is excluded. Even for reduced cross-sections, significant 
portions of the mass plane remain excluded. 

Acknowledgements 

We are greatly indebted to our technical collaborators, to the members of the CERN- 
SL Division for the excellent performance of the LEP collider, and to the funding agencies 
for their support in building and operating the DELPHI detector. 
We acknowledge in particular the support of 

Austrian Federal Ministry of Education, Science and Culture, GZ 616.364/2-III/2a/98, 
FNRS-FWO, Flanders Institute to encourage scientific and technological research in the 
industry (IWT), Belgium, 

FINEP, CNPq, CAPES, FUJB and FAPERJ, Brazil, 
Czech Ministry of Industry and Trade, GA CR 202/99/1362, 
Commission of the European Communities (DG XII), 
Direction des Sciences de la Matiere, CEA, France, 

Bundesministerium fur Bildung, Wissenschaft, Forschung und Techno logie, Germany, 
General Secretariat for Research and Technology, Greece, 

National Science Foundation (NWO) and Foundation for Research on Matter (FOM), 
The Netherlands, 
Norwegian Research Council, 

State Committee for Scientific Research, Poland, SPUB-M/CERN/PO3/DZ296/2000, 

SPUB-M/CERN/PO3/DZ297/2000, 2P03B 104 19 and 2P03B 69 23(2002-2004) 

FCT - Fundacao para a Ciencia e Tecnologia, Portugal, 

Vedecka grantova agentura MS SR, Slovakia, Nr. 95/5195/134, 

Ministry of Science and Technology of the Republic of Slovenia, 

CICYT, Spain, AEN99-0950 and AEN99-0761, 

The Swedish Research Council, 

Particle Physics and Astronomy Research Council, UK, 
Department of Energy, USA, DE-FG02-01ER41155. 
EEC RTN contract HPRN-CT-00292-2002. 



References 

[1] W. Loinaz and J. D. Wells, Phys. Lett. B445 (1998) 178; 

M. Carena, S. Mrenna and C. E. M. Wagner, Phys. Rev. D60 (1999) 075010; 

M. Carena, S. Heinemeyer, C. E. M. Wagner and G. Weiglein, hep-ph/9912223; 

M. Carena, J. Ellis, A. Pilaftsis and C. E. M. Wagner, Nucl. Phys. B586 (2000) 92. 
[2] J. F. Gunion, H. E Haber, G. L. Kane and S. Dawson, "The Higgs Hunter's Guide", 

Addison- Wesley Publishing Company, 1990. 
[3] X. Calmet and H. Fritzsch, Phys. Lett. B496 (2000) 190. 
[4] ALEPH Collaboration, A. Heister et al, Phys. Lett. B544 (2002) 25. 
[5] L3 Collaboration, P. Achard et a/.,Phys. Lett. B583 (2004) 14. 
[6] OPAL Collaboration, G. Abbiendi et a/., Phys. Lett. B597 (2004) 11; 

OPAL Collaboration, G. Abbiendi et a/., Eur. Phys. J. C40 (2005) 317. 



23 



[7] A. L. Read, in CERN Report 2000-005, p. 81 (2000). 

[8] E. Gross, G. Wolf and B. Kniehl, Z. Phys. C63 (1994) 417; 

Erratum Z. Phys. C66 (1995) 321. 
[9] P. Janot, in CERN Report 1996-01, Vol. 2, p. 309. 
[10] DELPHI Collaboration, P. Aarnio et al, Nucl. Instr. and Meth. A303 (1991) 233; 

DELPHI Collaboration, P. Abreu et al, Nucl. Instr. and Meth. A378 (1996) 57. 
[11] P. Chochula et al (DELPHI Silicon Tracker Group), Nucl. Instr. and Meth. A412 
(1998) 304. 

[12] S. Jadach, B. F. L. Ward and Z. Was, Comp. Phys. Comm. 130 (2000) 260. 
[13] E. Accomando and A. Ballestrero, Comp. Phys. Comm. 1999 (1997) 270; 

E. Accomando, A. Ballestrero and E. Maina, Comp. Phys. Comm. 150 (2003) 166. 
[14] A. Ballestrero et al, Comp. Phys. Comm. 152 (2003) 175. 
[15] T. Sj6strand et al, Comp. Phys. Comm. 135 (2001) 238. 
[16] F. A. Berends, R. Pittau and R. Kleiss, Comp. Phys. Comm. 85 (1995) 437. 
[17] P. Abreu et al, Nucl. Instr. and Meth. A427 (1999) 487. 
[18] S. Catani et al, Phys. Lett. B269 (1991) 432; 

N. Brown, W. J. Stirling, Z. Phys. C53 (1992) 629. 
[19] DELPHI Collaboration, J. Abdallah et al, Eur. Phys. J. C30 (2003) 447. 
[20] DELPHI Collaboration, J. Abdallah et al, 

Z~f* production in e + e~ interactions at y/s = 183-209 GeV 

(To be submitted to Eur. Phys. J. C). 
[21] DELPHI Collaboration, J. Abdallah et al, Eur. Phys. J. C32 (2004) 185. 
[22] A. Ballestrero et al, CERN Report 2000-009 (2000), p. 137. 
[23] DELPHI Collaboration, J. Abdallah et al, Eur. Phys. J. C32 (2004) 475. 
[24] T. G. M. Malmgren, Comp. Phys. Comm. 106 (1997) 230; 

T. G. M. Malmgren and K. E. Johansson, Nucl. Instr. and Meth. A403 (1998) 481. 
[25] G. C. Fox and S. Wolfram, Nucl. Phys. B149 (1979) 413. 
[26] T. Sj6strand, Comp. Phys. Comm. 28 (1983) 229. 

[27] DELPHI Collaboration, J. Abdallah et al, Eur. Phys. J. C32 (2004) 145. 



